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We study the anisotropic orbital effect in the electric field tunability of the phonon induced spin-flip rate 
in quantum dots (QDs). Our study shows that anisotropic gate potential enhances the spin- flip rate and 
reduces the level crossing point to a lower QDs radius due to the suppression of the Lande g-factor towards 
bulk crystal. In the range of 10'* — 10^ V/cm, the electric field tunability of the phonon induced spin-flip 
rate can be manipulated through strong Dresselhaus spin-orbit coupling. These results might assist the 
development of a spin based solid state quantum computer by manipulating spin-flip rate through spin-orbit 
coupling in a regime where the g-factor changes its sign. 



CD 

S 

a 

o 

> 
oo 
in 
in 
o 



X 



Controlling the single electron spins in QDs through 
the application of anisotropic gate potential is important 
for the design of solid state quantum computer Tun- 
ability of the phonon induced spin-flip rate and the elec- 
tron g-factor in III-V semiconductor QDs can be manip- 
ulated through the application of externally applied gate 
potentials4r— The strength of the Rashba-Dresselhaus 
spin-orbit coupling is determined by the asymmetric tri- 
angular quantum well along z-direction which approxi- 
mately estimates the density of electrons at the hetero- 
junction.'^'S Rashba spin-orbit coupling arises from the 
structural inversion asymmetry along the growth direc- 
tion and Dresselhaus spin-orbit coupling arises from bulk 
inversion asymmetry in the crystal latticci^ii^ 

Recently, it has been made possible to measure the 
electron spin states in gated QDs in presence of mag- 
netic fields along arbitrary direction.^ The physics be- 
hind this has been theoretically investigated by authors 
in Refs. [Ill and [l^ which confirms that the spin-orbit 
coupling can be used as a control parameter in the elec- 
tric field and magnetic field tunability of the electron g- 
factor tensor. The authors in Refs. [l3 and [13 measured 
the long spin relaxation times, approximately 0.85 ms in 
GaAs QDs by pulsed relaxation rate measurements and 
approximately 20 ms in InGaAs QDs by optical orienta- 
tion measurements. These spin-fiip rate measurements in 
QDs confirm the theoretical predictions of the suppres- 
sion of the phonon induced spin-flip rate by spin-orbit 
coupling with respect to the environment j^'*^'^'' Our work 
is along the lines of Refs. |3, [13, and [H with several new 
important findings. In particular, in this paper, we study 
the anisotropic orbital effect on the phonon induced spin- 
fiip rate for a system where the area of the symmetric and 
asymmetric QDs kept constant. Based on both theoreti- 
cal and finite element numerical simulation methods, we 
find that the anisotropic potential enhances the spin fiip 
rate and reduces the level crossing point to a lower quan- 
tum dot radius due to the suppression of the g-factor 
towards bulk crystal. 

We consider 2D anisotropic semiconductor QDs formed 
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FIG. 1. (Color online) Phonon induced spin-flip rate due to 
spin-orbit admixture mechanism as a function of QDs radius 
in symmetric QDs (a=b=l). We choose B=l Tesla. Also, 
inset plot shows the g-factor vs. QDs radius. The level cross- 
ing point occurs at = 73 nm. The material constants for 
GaAs QDs have been chosen from Refs. 13 and [l^ as follows: 
So = -0.44, m = 0.067, 7h = 4.4 A^ 70 = 26 eVA^, e/iw = 
2.34 X 10"^ erg/cm, = 5.14 x 10^ cm/s, st = 3.03 x 10^ cm/s 
and p = 5.3176 g/cm^. At £ = 7 x 10^ V/cm shown by 
dashed-dotted lines, the admixture mechanism due to spin- 
orbit coupling on the spin-flip rate is quite different because 
the electron spin states change their sign in these regime (see 
inset plot). 



in the conduction band in the presence of magnetic field B 
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along z-direction. The total Hamiltonian H — ±±xy- 
of an electron in the conduction band under the Kane 
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model^i^S, can be written as 

p2 1 fi 

Hxy = T, ^ ^"I'^ol"^;^ + by^) + -cr^w^, (1) 

2m 2 2 

-ffso = (o-x-Py - CTj/Pj;) + — {-a^Px + O-yPy) , (2) 



aR = 7fle£^, an = 0.787i 
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where P = — ifiV + is the 2D electron mo- 
mentum operator in the asymmetric gauge A — 
-^^-^{—yVb, Xyfa, 0) and uj^ — qq^lqB jh is the Zeeman 

frequency. Also, m is the effective mass, /is is the Bohr 
magneton, a is the Pauli spin matrices, ojo = is the 
strength of the parabolic confining potential with quan- 
tum dot radius ^o- The externally applied gate potential 
(V^y = l/2ma;Q(ax'^ -I- hy^)) in our theoretical model de- 
fines the lateral size of the QDs along x- and y-directions 
in the plane of 2DEG. By chosing = 10 — 100 nm, we 
mimic the experimentally reported lateral size of the QDs 
in Refs. [l] and[2l| and the potential induced in this range 
is much lower than the break down voltage (400 kV/cm) 
in GaAs heterojunctions. The strength of the Rashba 
and Dresselhaus spin-orbit couplings is determined by 
the relation an/aD = 1-5 x lO^^i?^/^ which tells us that 
ctR = at the electric field £' = 3 x 10^ V/cm. In the 
range of 10"* — 10^ V/cm, only the Dresselhaus spin-orbit 
coupling has an appreciable contribution in the manip- 
ulation of spin-hip rate in QDs. The asymmetric trian- 
gular quantum well potential {E — —dV/dz) arises along 
the growth direction and usually has a major contribu- 
tion to the E-filed tunability of the g-factor and spin- flip 
rate due to the interplay between Rashba-Dresselhaus 
spin-orbit couplings. The expression {2meE/h?)~^^^ es- 
timates the average thickness of the 2DEG where one can 
estimate the vertical average height of the QDs4.ii2, By 
chosing E — (0.1 — 10)10^ V/cm, we estimate the average 
height of the QDs from 2 nm to 10 nm which is in the 
range of experimentally reported values." 

The above Hamiltonian ([T]) can be exactly diagonal- 
ize d^^'^'^ and spin-orbit Hamiltonian can be used pertur- 
batively to find the energy states of the QDs. The energy 
spectrum can be written as 

Hxy — ("-+ + + 1) + (^4. — n_) fta;_ + — cr^w^, 

Hso = "i? (1 + i) [b^^'^K+ (s+ -i)a+ + h^^'^K+ (s_ -f i) a_ 
-|-a^/^?7_ {i - s_) a+ + a^/'^i^^ {i + s+) a_] 
(1 + i) [a^^^K- {i — s_) + a^^^n^ {i + S-|_) a_ 
-Ffoi/V (-^ + s+)a+ + b^/*i-i+ {i + s_) a_] + H.c, 
1 
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and Wc = ^ is the cyclotron frequency. 



We now turn to the calculation of the phonon induced 
spin relaxation rate in between two lowest energy states 
in QDs. The interaction between electron and piezo- 
phonon can be written aa^^i^ 



qa 



^i(q-r-a;<,ot)eAq„bqc 



2pVuJqa 



H.c. (4) 



Here, p is the crystal mass density, V is the volume of 
the QDs, creates an acoustic phonon with wave vec- 
tor q and polarization Bq, where a = l,ti,t2 are chosen 
as one longitudinal and two transverse modes of the in- 
duced phonon in the dots. Also, A^a = (jiQkePijke^^^ 
is the amplitude of the electric field created by phonon 
strain, where q = q/? and ePijk = ehu for i ^ 
^i* 7^ jij 7^ ^- The polarization directions of the 
induced phonon are e; = (sin 9 cos 0, sin 9 sin 0, cos 9) , 
gj J — (cos 9 cos 4>, cos 9 sin — sin 9) and = 
(— sin0, cos 0, 0). Based on the Fermi Golden Rule, the 
phonon induced spin transition rate in the QDs is given 
by4.5 
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27r 



/t^Ts E |M(q«)P^(fe.q-A), (5) 
J (27r) , , 



and H.c. represents the Hermitian conjugate. Also, 



where s;,St are the longitudinal and transverse acoustic 
phonon velocities in QDs. The matrix element M (qa) 
for the spin-flip rate between the Zeeman sublevels with 
the emission of phonon qa has been calculated pertur- 
batively.^"^ 

In Fig. m we quantify the influence of the Rashba- 
Dresselhaus spin-orbit admixture mechanism on the 
phonon induced spin-flip rate as a function of quantum 
dot radius for symmetric QDs. For the electric fields 
E = 10** and 10^ V/cm (shown by solid and dashed lines), 
we find that the transition time between spin up and 
down states increases with the increase in QDs radius and 
gives the level crossing point at £o = 73 nm. For the case 
E = 7x10^ V/cm (shown by dashed dotted line), the spin 
relaxation rate in QDs starts decreasing at £o — 35 nm 
because Zeeman spin splitting energy is very small which 
implies small phonon density of states. It becomes neg- 
ligible at £0 = 50 nm, and thus the spin relaxation rate 
turns to be zero. The spin relaxation time starts increas- 
ing at ^0 = 65 nm. However, in this regime, the opposite 
spin state is dominating because the g-factor of an elec- 
tron spin states changes their sign (see inset plots) . This 
is an important result for the design of spin based logic 
devices. Indeed, in GaAs/AlGaAs QDs, wavefunctions 
of electrons penetrate from GaAs QDs to AlGaAs bar- 
rier with the application of gate controlled electric fields 
where the g-factor of an electron changes its sigui^"— 
The level crossing takes into account, due to mixing, the 
Zeeman spin states |0, 0, — > and |0, 1, -I- > in QDs. The 
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FIG. 2. (Color online) Phonon induced spin-flip rate due to 
spin-orbit admixture mechanism as a function of QDs radius 
in asymmetric QDs (solid and dotted lines). As a reference, 
we also plotted spin-flip rate vs. QDs radius for symmetric 
QDs (dashed and dashed dotted lines). We choose the poten- 
tials characterized bya = 0.5&b = 2 for asymmetric QDs 
and a=b=l for symmetric QDs. Also we choose B = 1 T. 
As we see, spin-flip rate increases approximately by one half 
order of magnitude in asymmetric QDs. 
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FIG. 3. (Color online) (a) The anisotropic effect on the g- 
factor vs. QDs radius at the potentials characterized by a = 
b = 1 (solid line) for isotropic QDs and a = 0.5, 6 = 2 (dashed- 
dotted line) for anisotropic QDs. We choose E — 10^ V/cm 
and B = IT. Anisotropic potential gives the suppression of 
the g-factor towards bulk crystal and hence reduces the level 
crossing point to lower QDs radius. Accidential degeneracy 
appears in the range of 70 — 80 nm QDs radius which gives the 
cusp like structure in the spin-flip rate (see Refs. [TtI. IT^. and 
[25I ). (b) The interplay between Rashba and Dresselhaus spin- 
orbit couplings on the g-factor vs. the electric field in QDs 
induces the anisotropic effect due to the suppression of the 
g-factor towards bulk crystal. Here, we choose £0 = 20 nm, 
B = 1 T and a = 6 = 1. 



crossing point is theoretically investigated by the condi- 
^jQj^i7,i8 g^^^ _ ^ eo^i^+ i.e., h{uj+ - w_) = I^oImbS (see 
Eq. HJ. Substituting B = 1 T in the above condition, 
gives the crossing point at £q = 73 nm. Theoretically 
investigated level crossing point is in agreement with the 
numerically investigated values in the spin flip rate (see 
Fig. [1]). It can be seen that enhancement in the spin- 
flip rate occurs with the increase in electric fields. The 
level crossing point in the spin-flip rate is not aff'ected 
by the electric flelds which tells us that the level crossing 
point found in the spin flip rate is a purely orbital effect 
and is independent of the Rashba-Dresselhaus spin-orbit 
interaction. 

Fig. [5] explores the influence of anisotropic effects on 
the spin-flip rate vs. QD radius for the electric fields 
E = 10^ and E = 7 x 10^ V/cm. It can be seen that the 
anisotropic potential (a — 0.5, b — 2) enhances the spin- 
flip rate by approximately one half order of magnitude 
compared to that of symmetric potentials (a — b = 1). 
Note that we chose the above confining potentials in 
such a way that the area of the symmetric and asym- 
metric QDs are held constant. The level crossing point 
determined by the condition /i(w_|_ — w_) — jgolMs-B for 
anisotropic QDs is smaller than for the case of isotropic 
QDs if we held the area of the QDs constant. The cross- 
ing point for symmetric QDs was first studied by the 
authors of Refs. [l3 and llSl However, in this paper, we 
present the condition of the level crossing point for asym- 
metric QDs and by utilizing both theoretical and numer- 
ical methods, we report that the anisotropic potential 
reduces the level crossing point to a smaller QD radius 
as well as to smaller magnetic fields. Similar to Fig.[Tl at 
the electric field E — 7 x 10^ V/cm for anisotropic QDs 
in Fig. [2] (dotted lines), small Zeeman energy implies neg- 
ligible phonon density of states which gives zero spin fiip 
rate at £q = 46 nm and the g-factor with opposite sign 
(spin states change their sign) is observed at £0 = 63 
nm. Note that these numerically estimated values (zero 
spin-flip rate and the g-factor with opposite sign) occur 
at smaller QD radii for anisotropic QDs than in isotropic 
QDs. This tells us that the anisotropic potential leads to 
the quenching effect in the orbital angular momentunt^ 
that pushes the g-factor of an electron towards the bulk 
crystal which causes the level crossing point to occur at 
a smaller QD radius for anisotropic QDs (see Fig. [3ja)). 
From Fig.[31[b), we see that the g-factor can be tuned with 
the spin-orbit interactions. Recall that the electric held 
controls the strength of Rashba and Dresselhaus spin- 
orbit couplings (see Eq. [3]). Here, we again see that the 
suppression of g-factor towards bulk crystal induces the 
anisotropic effect due to the interplay between Rashba- 
Dresselhaus spin-orbit couplingSi ^^''^" 

To conclude, we have shown that the electron spin 
states in the phonon induced spin-flip rate can be ma- 
nipulated with the application of externally applied 
anisotropic gate potentials in QDs. The anisotropic po- 
tential causes the suppression of the g-factor towards bulk 
crystal that causes the enhancement of the spin-flip rate 



4 



and reduces the level crossing point to a smaller QD ra- 
dius. At sufficiently large electric fields, the phonon in- 
duced spin-flip rate can be tuned with Dresselhaus spin- 
orbit coupling by controlling the electron spin states in a 
regime where the g-factor changes its sign. 
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